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ABSTRACT 

It is now well established that stellar winds of hot stars are fragmentary and that the X-ray 
emission from stellar winds has a strong contribution from shocks in winds. Chandra high spectral 
resolution observations of line profiles of O and B stars have shown numerous properties that had 
not been expected. Here we suggest explanations by considering the X-rays as arising from bow 
shocks that occur where the stellar wind impacts on spherical clumps in the winds. We use an 
accurate and stable numerical hydrodynamical code to obtain steady-state physical conditions 
for the temperature and density structure in a bow shock. We use these solutions plus analytic 
approximations to interpret some major X-ray features: the simple power-law distribution of 
the observed emission measure derived from many hot star X-ray spectra and the wide range of 
ionization stages that appear to be present in X-ray sources throughout the winds. Also associated 
with the adiabatic cooling of the gas around a clump is a significant transverse velocity for the 
hot plasma flow around the clumps, and this can help to understand anomalies associated with 
observed line widths, and the differences in widths seen in stars with high and low mass-loss rates. 
The differences between bow shocks and the planar shocks that are often used for hot stars are 
discussed. We introduce an "on the shock" (OTSh) approximation that is useful for interpreting 
the X-rays and the consequences of clumps in hot star winds and elsewhere in astronomy. 

Subject headings: stars: early-type - stars: X-rays - stars: winds, outflows - stars: shocks - X-rays: 
stars 
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1. Introduction 

Our goal is to explain with clump bow shocks 
many of the problems found in the survey of 17 
normal OB stars by Waldron & Cassinelli (2007) 
(hereafter WC07) . Well- resolved spectral line pro- 
files are found to be neither shifted nor skewed to 
the blue (i.e., shortward) as had been expected 
(e.g., MacFarlane et al. 1991). The line widths 
are broader for supergiants and stars with thick 
winds than for lower luminosity stars. It seems 
clear that X-rays are formed in numerous shock 
fragments distributed throughout the wind. 

The topic of high non-radiative equilibrium 
temperatures in the winds of hot stars began with 
the discovery of superionization stages seen in the 
FUV spectra obtained with the Copernicus Satel- 
lite. Lamers & Morton (1976) analyzed the spec- 
trum of C Pup (04If) and found strong lines of 
O VI (1040 A) and N V (1240 A). These ioniza- 
tion stages seemed to require that the winds are 
"warm" with temperatures of about 2 x 10^ K. 
Alternatively, the ions could be produced by the 
Auger ionization process whereby 2 electrons are 
removed from the dominant ionization stage fol- 
lowing K-shell ionization by X-rays (Cassinelli, 
Lamers, & Castor 1978). Cassinelli & Olson 
(1979) used a "thin corona plus cool wind" model 
to show that the X-rays from a spatially thin 
corona plus the Auger process could explain ob- 
served boundaries in the HR diagram: a sharp 
cutoff in the presence of the O VI UV line that 
occurs at spectral type BO, a cutoff of the supe- 
rionization line of N V at B1.5, and similar ones 
for C IV at B5 and Si iv at B8 (Cassinelli & Ab- 
bott 1982). In each case the high ions could be 
explained by the removal of 2 electrons from the 
dominant stage of ionization, which are 0+'^ , 7V+^ , 
and C+^. Odegard & Cassinelli (1982) explained 
the more complicated case of the Si IV lines again 
using the Auger process. Soon after the predic- 
tions based on the Auger mechanism were made, 
X-rays were discovered from O stars in the first 
observations with the Einstein satellite (Seward 
et al. 1979; Harnden et al. 1979). However the 
observed X-ray spectral distribution did not agree 
with the idea that the X-rays were arising from a 
thin corona at the base of the cool wind, because 
the attenuation of soft X-rays by the cool wind 
was absent. It became clear that X-ray sources 



needed to occur farther out in the wind. Lucy 
& White (1982) and Lucy (1982) proposed that 
the X-rays were generated from shocks embedded 
in the wind and developed models of structured 
shocks. In fact, Lucy & White (1982) proposed 
that the shocks were at the outer face of clumps 
being driven through the wind, and thus were bow 
shock in nature. An advantage of the shock mod- 
els is that they could form naturally by virtue of 
the instability in the physics of line-driven winds 
(as first noted by Lucy & Solomon 1967). Thus 
unlike the warm wind and the corona plus cool 
wind models, shocks formed by line driven wind 
instabilities do not need to assume an input of me- 
chanical energy flux from the star to heat either 
an extensive part of the wind nor a thin coronal 
zone. 

Prom a moderately high resolution spectrum of 
Orion belt stars with the Einstein satellite Solid 
State Spectrometer, Cassinelli & Swank (1983) 
found that X-ray line emission from Si XIII and 
S XV in C, Ori. These indicated the presence of gas 
that is hotter than had been needed to explain the 
softer X-ray flux. Since this relatively hard line ra- 
diation could escape from deep in the wind, they 
proposed that there could still be zones, such as 
magnetic loops of very hot plasma, near the star's 
surface. Because single OB stars are not noted 
for X-ray variability at more than about the one 
percent level, Cassinelli & Swank also concluded 
that the embedded wind shocks could not be in the 
form of spherical shells, as in the picture of Lucy 
(1982), but rather there must be of order 10^ dis- 
tributed sources or "shock fragments" in the wind. 
The emission from a number of sources at a range 
of heights in the wind would lead to a statistically 
steady rate of X-ray production. Thus, even from 
relatively early in the history of hot star X-ray as- 
tronomy it was thought that the X-rays must arise 
from discrete distributed source regions, with the 
possibility that hotter zones are present near the 
star. Further improvements on that picture re- 
quired the higher spectral resolution of Chandra 
and XMM-Newton. 

The nature of the line-driven wind instability 
that leads to shocked regions was more fully ex- 
plained by Owocki, Castor, & Rybicki (1988), who 
developed a radial, 1-D picture for the spatial dis- 
tribution of the shocks. Spherically symmetric 
shock models have also been computed, for ex- 
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ample by Cooper & Owocki (1990) and Feldmeier 
et al. (1997). MacFarlane & Cassinelli (1989) de- 
veloped a basic model for an individual shock as 
being in the form of a "drivcn-wavc" with an in- 
ward boundary on the star-ward side of the wave 
where the wind collides with the driven wave, and 
an outward facing shock at the upper boundary 
where the driven wave catches up with slower mov- 
ing material ahead. The shock properties were ex- 
plained in terms of basic Rankinc-Hugoniot shock 
relations. The driven wave has a nearly constant 
pressure, set by the dynamic pressure of the in- 
cident material (pT^ci^), and the temperature of 
the gas can be quite high and is determined by 
Vrei^, where V^ei is the speed of the incident mate- 
rial relative to the shock front. For strong shocks 
the density increases by a factor of four relative 
to the incident material, and then in the layers 
of the drivcn-wavc where the gas radiativcly cools 
back to the radiative equilibrium value of the am- 
bient wind, the density increases by about 3 or- 
ders of magnitude. It is likely that a shell of such 
a high density contrast gas is unstable, such that 
the compressed gas will not remain in the form of 
a spherically expanding shell, but instead break 
up into clumps with densities of about 10'^ times 
that of the ambient wind. 

The initiation of our bow shock approach arose 
through our attempt to understand the anomalous 
properties of the BO. 5 V star r Sco (Howk et al. 
2000). In Copernicus spectra (Lamers & Roger- 
son, 1978), the star shows unusual red-shifted ab- 
sorption features in the O vi and N v P-Cygni 
lines. This could arise if there were infall of mat- 
ter. In addition, ROSAT observations of r Sco 
showed that the X-ray spectrum is unusually hard 
relative to other early B stars. Howk et al. pro- 
posed that clumps form in the wind, become dy- 
namically un-coupled from the flow and line forces, 
and if they have a certain range in clump mass, 
follow trajectories whereby they fall back toward 
the star. The bow shock around the clump would 
account for the X-ray emission, and drag owing to 
the wind/clump velocity difference would be im- 
portant in the clump trajectory. 

The idea that hot star winds are clumpy has 
been especially well investigated in regards to 
Wolf-Rayet (WR) stars. Lepine & Moffat (1999) 
found that emission features seen moving across 
broad optical lines of He II lines can be ex- 



plained as clump emission regions moving out in 
the winds. Hillier (1991) showed that unclumped 
winds of WR stars would lead to line broadening 
by electron scattering that is not observed. Nugis 
& Lamers (2000) found that the observed IR flux 
distribution of WR stars is better flt with clumped 
models. There is also evidence for clumping in the 
case of other early-type stars. Lupie & Nordsieck 
(1980) found that the position angle of polariza- 
tion seen in B supergiants shows irregular changes, 
and they interpreted this as arising from very large 
density enhancements in the winds. This was per- 
haps the first indication of fragmentation at a large 
scale. Brown, Ignace, & Cassinelli (2000) showed 
that only large scale concentrations of matter can 
lead to observable polarization changes, because 
smaller scale reshuffling of electrons in an enve- 
lope does not lead to a significant net polariza- 
tion change. In addition to the WR and B su- 
pergiants, there is new evidence of wind clumping 
in O star winds (Bouret et al. 2003; Evans et al. 
2004; Bouret, Lanz, & HiUier 2005). 

Clumping is important a) because it leads to 
overestimates of mass-loss rate (M) values of 
early-type stars that are derived from density 
square diagnostics (Abbott, Bieging, & Church- 
well 1981), and b) because the downwardly re- 
vised M values affect the X-ray spectra owing 
to a reduced absorption column density through 
overlying wind material. The best estimates of M 
had long been assumed to be those derived from 
the free-free radio flux of stars (Barlow & Cohen 
1976; CassinelU & Hartmann 1977). The results 
were considered the most reliable because the free- 
free transitions lead to a large (LTE) opacity; the 
observable radio flux tends to form in the outer 
regions of the wind where the velocity is constant 
and the density p varies simply as r~^; and the 
wind temperature cancels from the emergent flux 
formula. 

The topic of clumping has recently become a 
major one because surprisingly large reductions to 
the mass-loss rates for O stars have been suggested 
by Massa et al. (2003) and FuUerton, Massa, & 
Prinja (2006). These authors argue that a bet- 
ter estimator of M is the P-Cygni profile of the 
P V (All 18, 1128A) doublet. This spectral feature 
has several advantages: it originates in a dominant 
ion stage, so there is a minimal ionization fraction 
correction and the profile is typically unsaturated 
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in UV spectra. Also, since the line opacity de- 
pends only linearly on density, the line depth is 
not affected by clumping in the wind. Fullerton 
et al. deduce mass-loss rate reductions by at least 
a factor of 1/10 to 1/100 times traditional values. 
However, based on our clump picture, we ques- 
tion the idea that the clumping does not affect 
the abundance of the P V, because the clumps are 
immediately adjacent to X-ray sources, and thus 
the fraction of phosphorus in the observable P V 
stage could be reduced significantly owing to the 
multiple ionizations associated with the Auger ef- 
fect (Odegard & Cassinelli 1982) 

Another argument against the significant de- 
crease in M is from WC07 who find that the radii 
of formation of the X-ray lines corresponds quite 
well with the values at which radial optical depth 
is unity when using the traditional M values. Thus 
lines from soft ions such as that of N vil (25 A) 
form at 5 to 10 stellar radii because the opacity to 
radiation at these wavelengths is large. (Contin- 
uum opacity, much of it K-shell opacity of abun- 
dant metals, varies roughly as A^.) Shorter wave- 
length lines, that originate from the higher ions, 
such as Ne IX, Mg xi, Si xiii, can be seen as form- 
ing progressively deeper. These are also near their 
respective optical depth unity locations, if tradi- 
tional M values are used. Another good argument 
against the drastic decrease in the mass-loss rates 
are that it would aflfect the well-established results 
of massive star evolution (Hirschi 2007). The topic 
of clumping is an active one and at the recent 
Potsdam international workshop on hot star wind 
clumps (Hamann, Feldmeier, & Oskinova 2007), 
a consensus was reached that M could not be re- 
duced by more than about a factor of 3 from the 
traditional values. 

In summary, it has become important to un- 
derstand the properties and the effects of clumps 
in winds. Our plan here is to consider the various 
effects discrete clumps would have on observations 
of X-ray lines and to establish an analytical tool 
to interpret observations. There is evidence for 
clumps being important in astrophysics in general, 
for example, bow shocks also appear in images of 
Herbig-Haro objects (e.g., Eisloffel et al. 1994) and 
planetary nebulae (Odell et al. 1995). By using 
Chandra data we can obtain good information on 
the clumps in hot star winds. 

In §2 we summarize some of the problematic 



X-ray results obtained from Chandra. In §3 we 
describe the numerical method we have used and 
show the most relevant results regarding the post 
shock temperature and density structure as well 
as the temperature distributed emission measure. 
In §4, the bow shock structure is contrasted with 
plane parallel shock pictures and in §5, we develop 
a simple "on the shock" approximation to provide 
some insight regarding bow shocks. In particular 
wc give a derivation of the distributed emission 
measure power law with temperature and discuss 
the angular flow around the clump. Overall con- 
clusions are summarized in §6. 

2. Hot Star X-ray Problems 

At the forefront of current problems in this field 
is the finding that X-ray line profiles from massive 
star winds are quite different from what had been 
expected. MacFarlane et al. (1991) had predicted 
that lines formed by shocks in stellar winds should 
be blue shifted and skewed. This is because the 
shock regions on the far side of the star (where 
the gas is red shifted) would be more highly at- 
tenuated by bound-free continuum opacity of the 
wind matter by virtue of being at higher column 
density relative to the observer. This effect was 
below the spectral resolution of the emission lines 
of the two B stars (e CMa and /3 CMa) that were 
observable with the Extreme Ultraviolet Explorer 
{EUVE) satellite. These stars were observable at 
70-730 A and 500-700 A, respectively (CassinelU 
et al. 1995, 1996). However, well-resolved line pro- 
files at the even shorter X-ray wavelengths are 
observable with the high energy and medium en- 
ergy grating spectrometers (HEG and MEG) on 
Chandra. The Chandra observations (Waldron 
fc Cassinelli 2001; Miller et al. 2003; Leuteneg- 
ger et al. 2006; WC07) do not show the expected 
blue-ward skewing of the lines nor blue shifting 
of the line centroids. The lines tend also to be 
broad, although much narrower for the low lumi- 
nosity classes of hot stars, and for all luminosity 
classes the HWHM is less than the terminal wind 
speed. This problem of symmetry alone has moti- 
vated several ideas. 

A) If the mass-loss rates were reduced by an order 
of magnitude or more, the winds would be 
suflaciently thin that radiation from both the 
near and far side could escape (Cohen et al. 



4 



2006; Leutenegger et al. 2006). 

B) Owing to clumps, the winds could be more 

porous to the transfer of X-rays (Oskinova 
et al. 2006). The presence of high density 
fragments allows the X-rays to escape from 
deeper in the wind, and one could detect suf- 
ficient X-rays from the far side of the star 
to eliminate the net shifting and skewing ef- 
fects. 

C) There is enhanced Sobolev escape of line radi- 

ation out the side of shocks (Ignace & Gayley 
2002). This means that the line radiation is 
coming in large part from the "sides" of the 
star, as seen by the observer and this occurs 
where the radial flow has a small line-of-sight 
velocity. Although we will not be employing 
the Sobolev escape argument, we will find 
the enhanced sideways escape to be a useful 
concept. 

In the study of clump effects by Feldmeier et al. 
(2003) and Oskinova et al. (2004), the clump re- 
gions are pictured to be in the form of plane par- 
allel slabs. These correspond to regions of en- 
hanced absorption of the X-rays that originate in 
broader more diffuse wind regions. WC07 note 
that the derived radii of the line source regions 
found by these authors are larger than the rather 
small radii that WC07 find from an analysis of the 
forbidden-intercombination-resonance (fir) triplet 
lines of He-like ions. Oskinova et al. (2004) have 
developed plausible explanations for the unshifted 
un-skewcd lines that are observed in C Orionis. 
The clumps in their view arise from a runaway 
effect first described in Feldmeier et al. (1998). 
Whether porosity at the required level can explain 
the observations has been questioned by Owocki & 
Cohen (2007). Nevertheless, the Oskinova et al. 
picture combines several elements: the fragmen- 
tary nature of shocks and the possibility of pro- 
ducing X-rays at a range of radii that is needed in 
any explanation of hot star X-rays. 

3. A Hydrodynamical Calculation of a 
Wind Colliding with an Impenetrable 
Object 

We use the magnetohydrodynamical (MHD) 
code that was developed by Cho & Lazarian (2002) 



to model the bow shock that forms from a plane- 
parallel flow impinging upon a spherical blunt 
obstacle. As a flrst step in understanding the 
effects of bow shocks on hot star X-ray emis- 
sion, we choose to assume that the clump is 
impenetrable. Actual clumps are likely to be 
more complicated, and the clumps coiild have a 
highly time-dependent interaction with the inci- 
dent wind. However to ensure that we are not 
dealing with numerical noise phenomena, we treat 
the simplest case possible as a starting point for 
our theoretical exploration of bow shock effects. 
The main expectation is the formation of a bow 
shock around the clump that exhibits a range in 
hot plasma temperatures along with a non-trivial 
vector velocity flow. This recipe is suggestively 
promising for explaining the observed X-ray spec- 
tral features from hot stars. 

For the simulation, the clump is assumed sta- 
tionary and spherical with a radius Re- The inci- 
dent flow is plane-parallel at constant speed and is 
parameterized by the mass flux and flow speed. In 
such an idealized case, the X-ray emission arises 
owing to the post-shock gas that envelopes the 
clump. All of the emission arises from the wind 
matter that is heated, compressed, and redirected 
at the bow shock interface. For the purpose of 
interpreting X-ray observations, we are interested 
in flnding the temperature and density structure, 
and the distributed emission measure. Using these 
physical properties in conjunction with a cooling 
function provides the X-ray emissivity of a blob 
as a function of wavelength and the line proflle 
emission for a particular transition of interest. 

For computing the structTire from bow shocks, 
we use a third-order hybrid essentially non- 
oscillatory (ENO) MHD routine as described in 
Cho & Lazarian (2002) and references therein 
to reduce spurious oscillations near shocks, two 
ENO schemes are combined. When variables are 
sufficiently smooth, the 3rd-order Weighted ENO 
scheme (Jiang & Wu 1999) is used without charac- 
teristic mode decomposition. When the opposite 
is true, the 3rd-order Convex ENO scheme (Liu & 
Osher 1998) is used. A three-stage Runge-Kutta 
method is employed for time integration. The 
ideal MHD equations are solved with magnetic 
field set to zero everywhere. The simulation runs 
until it reaches an approximately steady-state bow 
shock, usually after about two wind crossing times 
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owing to the highly hypersonic nature of the flow. 
The simulations reported in this paper are axisym- 
metric. We adopt cylindrical coordinates (ro, (jj, z) 
with incident flow moving in the +2;-direction and 
so the bow shock is axisymmetric about the z-axis. 
We shall frequently refer to the transverse coor- 
dinate w as the impact parameter. Some early 
results were presented by Moeckel et al. (2002). 

In the calculations a polytropic relation with 
P (X p''' is assumed so that the temperature struc- 
ture is determined by the expansion cooling of the 
gas after it passes through the bow shock geome- 
try. Adiabatic cooling is an appropriate approxi- 
mation for the case of high velocity inflow, so the 
post-shock region is at a high temperature where 
radiative losses are less important than expansion. 
The computed results for the shape of the shock 
are quite similar to the semi-analytic results of Lo- 
max & Inouye (1963) for a 7 = 5/3 bow shock. 

The incident velocity is specified as a Mach 
Number Mw, and here we present the results of 
an analysis of an adiabatic bow shock with wind 
flow Mw = 47 and 71, that correspond to relative 
velocities of 1000 and 1500 kms^^, respectively. 
Throughout this paper, the subscript 'W refers 
to wind properties. A fast wind is assumed to be 
incident upon a rigid stationary sphere, creating 
the bow shock structure. The sphere has a ra- 
dius of 32 grid-points. The center of the sphere 
was located one half grid spacing below the lowest 
grid-point. The results of this simulation are ro- 
bust for other large values of Mw since the shock 
shape, flow pattern, and shock structure are all 
nearly independent of Mach number in the hy- 
personic limit of Mw ^ 1 for an adiabatic case 
(Hayes & Probstein, 1966). The parameters used 
in our simulations are given in Table [TJ 

The most basic parameters for determining the 
structure of a bow shock around a blunt object 
are the incident velocity, which determines the 
post-shock temperature, and the incident mass 
flux pv^K-ei- The numerical calculation is made 
in the rest frame of the rigid sphere where Kei 
is the difference between the radial wind veloc- 
ity and the radial clump velocity. The post-shock 
gas temperature follows the well known Rankine- 
Hugoniot relation between the velocity perpendic- 
ular to the shock and the post-shock temperature. 
The maximum temperature, Ta achieved occurs 
at the "apex" of the bow shock (i.e., along the 



line of symmetry) and is given by 

= UMk( ^^'^ , V (2) 

Viooo kms-iy ^ ' 

where in the latter expression we have evaluated 
the constants using fully ionized solar abundances 
and mu — 0.62, and VV,_l is the perpendicular 
speed of the pre-shock gas relative to the shock 
front. 

Figure[T]shows simulation results from our adia- 
batic bow shock case using the parameters listed in 
Table [T] (using the second set of parameters wher- 
ever two values are listed). Figure [T] shows (a) 
streamlines of the flow, (b) the density, (c) tem- 
perature, and (d) emission measure distributions . 
Figure [2] shows temperature contours behind the 
bow shock superposed with vectors detailing the 
post-shock flow velocity. Note especially the de- 
creasing temperature along the bow shock with 
distance from the apex, and also that the flow of 
gas around the clump leads to a significant trans- 
verse vector velocity field in conjunction with the 
temperature distribution. Note also, that all the 
X-ray emission properties, in which we have spe- 
cial interest here, closely hug the bow shock sur- 
face. This is what inspires our introduction of the 
"on the shock" or OTSh approximation. In the 
following section, we consider observables relating 
to these properties of the clump- wind simulation. 

4. The Temperature Distribution of the 
Emission Measure 

4.1. Results for a Single Clump 

An important property for understanding pro- 
cesses that generate hot plasma in astrophysi- 
cal sources is the temperature distribution of the 
emission measure. This distribution is often called 
the "differential emission measure" , but the use of 
this term in the X-ray literature is not uniform, 
and we prefer simply to refer to the temperature 
distribution of the emission measure or EM{T). 
The contribution to the emission measure at a 
given temperature range from T to T + dT can 
arise from totally disconnected regions in the X- 
ray emitting region. The volume emission measure 
(EM) is defined as 
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EM = J N^NpdV, (3) 

a volume integral involving the product of the elec- 
tron and proton number densities iVc and Np. Wo- 
jdowski & Schulz (2004, 2005; hereafter WS04 and 
WS05) find that the amount of hot plasma at each 
range in temperature to be a decreasing power law 
for most hot stars, with the exception of those 
few hot stars known to possess highly magnetized 
envelopes (e.g., 0^ Ori C, see Donati et al. 2002). 
Each of the X-ray emission lines is associated with 
a range of temperatures that can be found from 
the APEC software (Smith et al. 2001). From the 
strength of the line, Wojdowski & Schulz find the 
interesting result that EM{T) is a downward slop- 
ing power law extending from about 2 MK, up 
to an apparent maximum that is typically about 
20 MK. Such a distribution holds for the majority 
of the stars that were analyzed. For stars known 
to be highly magnetic, the EAI{T) distribution 
can have a very different shape. The highly mag- 
netic star 0^ Ori C even has a positive slope over 
the whole temperature range indicating the pres- 
ence of larger amounts of increasingly hot mate- 
rial. Clearly such a distribution would have to 
stop at some high temperature, although this is 
not seen in the 0^ Ori C plots shown by WS05 or 
WS07. Table [2] lists values of the EM{T) power 
law index for several OB stars as derived from the 
analyses of WS05 and WC07. 

Using the hydrodynamical code, we obtain a 
power-law emission measure distribution EM{T) k, 
lO^^ "* cm-3 X {T/Ta)'P, shown in Figure [H where 
p is approximately —4/3. Note that the drop-off 
of at low temperatures is an artifact of the cal- 
culation that arises from the fact that the grid 
in the impact parameter direction is truncated at 
some Umax- When this maximal impact parame- 
ter is increased the slope continues to follow the 
—4/3 slope to lower temperature values. With 
this power-law distribution, the temperature gra- 
dient dEM/dT oc TP-i oc T^'^/^, a result we wih 
use later on. This figure was obtained by binning 
the emission measure results of each simulation 
volume element into temperature cells of 0.1 in 
logT. 

It is useful to have the following scaling of the 
results of the simulations to parameters for a dif- 
ferent clump: with radius R, relative wind/clump 



velocity Vw, wind density Nw {— y^NcNp), and 
apex temperature Ta. We derive the scaling 

dEM _ lOSi'^cm-^ 
dT ~ 14 MK ® 

where Vg is the incident speed in 10* cm/s. Here 
we have used subscript, s, for our simulation re- 
sults. The -p — —4/3 slope in Figure [3] somewhat 
underestimates the emission from the highest tem- 
perature bin in the distribution (by a factor of 
about 2) , and we discuss this later in our analytic 
derivation of the dEM/dT relation. 

From an interpretational perspective, the most 
important feature to note about the power law is 
that there is significantly more gas at low tempera- 
tures than at the high temperatures near the apex 
of the shock. The increasing emission measure 
for the lower temperatures (and hence the lower 
ion stages) arises from the fact that a bow shock 
has an increased area toward the wings where the 
shock is increasingly oblique. As a result of the 
power law distribution EAI{T), clumps that are 
deep in a wind can produce a great deal of low ion 
emission, say at O Vli. In the case of stars with 
high mass-loss rates, such as Of stars and OB su- 
pergiants, the wind opacity can substantially block 
this emission from direct observations. 

The stars that have been studied most inten- 
sively in the history of wind theory and hot star X- 
ray astronomy are C Ori and C Pup, and both have 
EM{T) slopes reasonably close to the p = —4/3 
value. However some stars, such as r Sco, have 
a very shallow but negative slope. It has recently 
been discovered that this particular star is also one 
with a strong magnetic field (Donati et al. 2006), 
thus it is likely that this star has a combination 
of X-rays from both the clump regions and the 
positive sloping EM{T) magnetic regions. Sort- 
ing the contributions from each could be difhcult; 
however, it appears that the EM{T) analyses can 
yield valuable new information about clump prop- 
erties and may also provide an indicator of fields 
in stars that have not yet shown measurable Zee- 
man effects. We have found a case of a positive 
sloping EM{T) in our recent calculations of X- 
rays from Be stars (Li et al. 2007) using a Mag- 
netically Torqued Disk model for Be stars. This 
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model is based on the idea that the matter that 
enters a Be star disk is from a wind that has been 
channeled by a dipolar magnetic field and which 
is co-rotating with the star (CassincUi ct al. 2002; 
Brown et al. 2004). Higher temperature gas origi- 
nates farther from the star since the impact of the 
wind with the disk is at a higher speed. In this 
example there is a maximum to the EM(T) dis- 
tribution that is directly related to the magnetic 
field at the base of the wind. Next we consider the 
consequences of having a large number of clumps 
contributing to the EM{T). 

4.2. Ensemble of Clumps 

An inspection of the results from WS05 reveals 
that, except for the magnetically dominated star 
6'^ Ori C, observed EM{T) distributions have neg- 
ative power-law slopes. However the values are not 
exactly the —4/3 that we derive from the clump 
simulation. This can arise for several reasons. 
The first one has already been discussed, namely 
some stars are strongly magnetic and those fields 
can dominate the EM (T) following physics that is 
quite distinct from clump bowshocks, and so the 
results are not surprisingly quite different. How- 
ever the star HD 206267A (O 6.5V(f)) has a slope 
that is steeper than the —4/3 value. In WC07 we 
found that for near main-sequence stars, one could 
detect source emission coming from a wide range 
in depths as the winds are not as optically thick 
as those in the supergiants. Thus a second way 
to explain a deviation from a —4/3 slope in the 
case of steeper EM{T) is to consider the X-rays 
as arising from an ensemble of different clumps. 

Such an ensemble docs not represent only more 
or fewer clump bowshocks but in fact a range of 
Ta values. Each individual clump contributes a 
EM{T) distribution of slope -4/3 for T < Ta, 
but now clumps exist at different radii and thus 
a range in Ta may plausibly exist from a range 
in Vic\ values. The consequence of staggered val- 
ues in Ta leads naturally to a slope that is steeper 
than —4/3. All clumps contribute to low tempera- 
ture plasma, but only a small minority contribute 
to the absolute maximum temperature achieved in 
the entire wind. Note that one way to provide an 
observational cut-off for lower temperatures is to 
recognize that even rather low mass-loss winds will 
eventually lead to substantial photoabsorption at 
sufiiciently low X-ray energies. Second, our simu- 



lation assumes the clump is small compared to the 
radius of the star, such that the wind is plane par- 
allel on the scale of the clump size. However, the 
stellar wind is in fact spherically divergent, and so 
the flow striking the a large bow shock is actually 
more oblique than is achieved in the plane-parallel 
flow simulation. 

4.3. Adiabatic and Radiative Cooling Re- 
gions 

Before leaving the topic of EM{T), it is impor- 
tant to comment on the applicability of the adia- 
batic assumption. Across a bow shock there is a 
large range in temperatures, and even if the adia- 
batic assumption holds near the apex of the shock, 
it will fail somewhere out in the wings, since radia- 
tive cooling is extremely efficient for T < lO^K. 
So it is relevant to consider the regimes in which 
the cooling time is small or large compared to the 
flow time in the simulation. 

Adiabatic cooling is a result of expansion cool- 
ing as the flow navigates around the clump. 
Throughout this trajectory the gas is emitting 
X-rays, and so the gas is cooling by radiation as 
well; however, the adiabatic assumption adopted 
in our hydrodynamical simulation will apply to 
those cases when the radiative cooling timescale 
is slow relative to expansion cooling. 

To characterize the comparison, we define ffiow 
as the flow time scale associated with adiabatic 
cooling. This will scale as tnow ~ Rc/Vroi- With 
Rc = 10^° cm and a wind speed of 500 kms~^ 
(approximately a quarter of the terminal speed 
for a typical O star wind), the flow time works 
out to <flow ~ 10^ s. We define the radiative 
cooling time with trad, which can be estimated as 
the ratio of the thermal energy density Uth to the 
cooling rate dU/dt. Thus, trad ~ Uth/{dU/dt) = 
{NekT)/[A{T)N^], where A is the cooling function 
in erg cm^ s~^ which is tabulated in Cox (1999), 
giving the result log A = —21.6 — 0.6 (logT — 5.5). 
At a value of T = 10^ K associated with Ta and 
a number density of A^e = 10^° cm~^, the radia- 
tive cooling time becomes trad ~ s, about ten 
times larger than the flow timescale under these 
conditions. 

Clearly, the radiative cooling is dominant over 

adiabatic cooling where the wind density is large 
or the shock temperature is low, and in future sim- 
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ulations it will be important to include radiative 
cooling. However, it is still useful to consider the 
the limiting case of adiabatic cooling. First, al- 
though ID simulations show radiative cooling to 
be important, expansion cooling is there limited 
to divergence, and shock structures are neces- 
sarily spherical shells in these simulations, which 
we know not to be true based on low levels of X- 
ray variability (Cassinelli & Swank 1983). Our 
2D simulations allow for a new geometric avenue 
of expansion cooling because of the channeling of 
the gas around the clump in the form of a bow- 
shock, thus requiring a new assessment of how adi- 
abatic cooling contributes to the interpretation of 
hot star X-rays. Second, in the hypersonic limit, 
the results for pure adiabatic cooling are quite ro- 
bust because the bowshock shape and tempera- 
ture distribution are independent of density. This 
is not the case for radiative cooling, hence one ex- 
pects the radiative bowshock properties to be dif- 
ferent for clumps at different radii in the wind. 
And third, as discussed in the preceding section, 
we can explore the extent to which this limiting 
case can reproduce observed differential emission 
measure distributions. Departures from the model 
predictions are suggestive of contributions by ra- 
diative cooling or other effects, such as stellar mag- 
netism. 



wind flow relative to the center of the clump. The 
shape of the bow shock from our simulation can 
be fit with a power-law curve given by 

with a — 0.35 and m = 2.34, hence a shape not 
far from a parabola. This provides a good fit 
out to w ^ clump radii. An exact solution to 
the shock shape derived using an inverse method 
is given in Lomax & Inouye (1963), and our re- 
sult for the shock shape agrees well with theirs. 
Their paper does not provide the temperature and 
emission measure information of special interest to 
us. In addition to the shape, it will also be use- 
ful to know the derivative of the shape (i.e., the 
position-dependent tangent). We define this to be 
g{'Dj) = t&uAi = dz/dw where Ai is the angle 
that the incident wind makes relative to the bow 
shock. 

5.1. Post-shock velocity components 

Figure |4] shows a post-shock (hereafter identi- 
fied with subscript "P" ) flow trajectory associated 
with crossing the bow shock. The jump condi- 
tions for the velocity components perpendicular 
and parallel to the shock front are: 



5. The "On the Shock" Approximation 

From the numerical modeling, we find, in con- 
trast with planar shocks (e.g., Feldmeier et al. 
1997), that the X-ray emission is dominated by a 
zone lying just behind the shock front (as seen in 
Fig. [Hi) . Hence, we introduce the "On the Shock" 
or OTSh approximation for a simple analysis of 
bow shocks. In the case of a planar frontal shocks, 
the post-shock flow must continue in a straight 
line, and thus the only cooling that occurs is ra- 
diative emission. In the case of a bow shock, there 
is always an expansion of the gas around the sides 
of the clump and is associated with adiabatic cool- 
ing. 

Consider a simple but useful picture for the 
shock structure. As before, let the z axis be the 
radial direction from the star through the apex of 
the shock and through the center of the spheri- 
cal clump. Let vj be the perpendicular direction 
that corresponds to the impact parameter of the 



Vp^i_ = -Vw.± (6) 

yp,\\ = yw,\\ (7) 

All across the face of the shock the incident wind 
speed is given by V = Vw^z- Based on the ge- 
ometry of Figure |4] we can derive the following 
relations: 

W._L = Vw cos(Ai) = Vw J- — - (8) 

VV,|| = Vw sin(Ai) = Vw - J== (9) 

where we have expressed the velocity quantities in 
terms of the shape gradient g{vj). The following 
limiting behavior is implied: as g <C 1, the wind 
is essentially normal to the shock as occurs along 
the stagnation line, and at large impact parameter 
where g ^ 1, the shock becomes nearly parallel to 
the incident wind flow. 
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Consider a jump that occurs at some position 
on the shock {w, z). The incident wind speed Vw 
relative to the orientation of the shock front can 



be expressed as Vw — Vw.x' 



The total 



post-shock velocity is Vp — (Vm/,±/4)^ + Vvi/,||' 
which upon using substition using the relations ([8 
and ini) becomes 



(10) 



[1 + /Z?] ' 

where on the post-shock side, the velocity vector 
makes an angle relative to the normal implicitly 
defined by h\ = tanA2 = Vp^±/Vp^\\ = 4 tanAi. 
Thus, 



Vp__L = Vp cosyl2 = Vp 
Vp II = Vp sin A2 — Vp 



hi 



(11) 

(12) 



5.2. Temperature distribution across the 
bow shock 

The post-shock temperature as a function of 
impact parameter w is related to the change in the 
perpendicular velocity component as described in 
the previous section. The temperature distribu- 
tion is given by 



rp t \ ^ ^"^^ (\T ' ^2 

Tp{w) = 16~A~^ "-^ 
= Ta cos2(Ai) 

l+g\uj) 



(13) 
(14) 



where Ta is the highest temperature gas, at the 
bow shock apex. Thus we see that each impact 
parameter point on the shock has a specific post- 
shock temperature associated with it, hence for- 
mally vj = f{T). The simple power law for the 
bow shock geometry in ([5]) leads to 



— =g{vj)=am{m/R,y^-\ 



(15) 



thus, we get an explicit equation for temperature 
in the OTSh approach: 



rp / \ 2.68 

[!+/!-'= l + 0.67(£ 

(16) 

This expression is plotted against the maximum 
temperature along an impact parameter in Fig- 
ure [5] where the agreement between the model and 
the preceding expression is remarkably good. This 
close agreement is related to the conditions of the 
simulation being hypersonic. 

5.3. Polar angle distributions of the X- 
rays 

To facilitate the use of these results for eval- 
uating observables from an ensemble of clumps, 
such as in predicting emergent line widths, it is 
convenient to express the angles of the inflow and 
post-shock flow relative to the local "radial" (or 
z) direction. This is a blob-centered (r, 9, (p) sys- 
tem. The incident radial wind flow has 9in — 0°. 
On the post-shock side, the angle is 6p = A2 — Ai, 
hence the stream flow emerges from the bow shock 
at 



tanOp 



3.9 



(17) 



1 + 452 - " 

Thus the "transverse" or sine component of the 
post-shock velocity vector is sin = /12/ ■\/l + 

What is particularly interesting about the post- 
shock velocity and its direction is that it affects the 
observed line profiles. The post-shock speed can, 
for example, have a rather large transverse value, 
as illustrated in Figure [H To choose one case for 
example, consider the one angle quadrature value 
of Lucy & White (1981): Ai = 30° then A2 = 67°, 
Vp = 0.56Vw and the transverse component of the 
post-shock velocity is one-third of Vw- This is al- 
ready a significant fraction of observed X-ray line 
widths (in WC07 this was typically about 0.4Vjv). 
In a simulation involving many clumps, one would 
need to account for the line-of-sight velocity of 
each of the clumps. The narrower lines found by 
WC07 for the near main sequence stars and for 
the lower ion stages can be explained by having 
the line formation regions being closer to the star 
where the local wind speed is smaller than is the 
case for the supergiants. For the latter case the 
radiation must arise from regions above about op- 
tical depth unity, and this is well out in the wind 
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where the wind speed is near the fast terminal 
speed. 

The transverse component of the X-ray emit- 
ting matter is another significant difference from 
planar shocks. The planar shocks are usually as- 
sumed to be perpendicular to the flow, and thus 
the post-shock matter continues flowing only per- 
pendicular to the shock front. 

5.4. Temperature distributions of the emis- 
sion measure 

Semi-analytic expressions for bow shock prop- 
erties have been derived before, especially for 
the simple case that we consider here in which 
the compact object is an impenetrable sphere. 
For example Canto & Raga (1998) derived the 
shock shape and pressure density structure for an 
isothermal bow shock. This is valid for cases in 
which the cooling time is rapid and so the shock 
material forms a thin sheet around the blunt ob- 
ject. In contrast our specific interest is in the hot 
adiabatic region where X-ray emission originates. 

Our results for the post-shock density have mo- 
tivated our "on the shock" approximation that 
provides a temperature, velocity, and EM from 
each spatial element of the bow shock surface. 
These are useful for understanding some of the 
outstanding questions raised by high spectral res- 
olution X-ray data (e.g.. Burke et al. 2006). In 
a separate paper, the OTSh approximation will 
be employed in quantitative spectral synthesis to 
model X-ray line profile shapes. Here, we present 
a derivation for EM{T) to show that the power- 
law temperature distribution arises naturally from 
the bow shock geometry. 

Based on our hydrodynamic simulation, the 
OTSh approximation asserts that the hot gas con- 
tributing to EM closely hugs the bow shock sur- 
face. In this region the density is nearly constant 
at = 'iNw everywhere along the bow shock. 
Hence, the volume element, dV , associated with 
dEM must scale roughly as the incremental sur- 
face area dS of the bow shock times some depth 
parameter which we call A£. This A£ parameter 
represents a characteristic depth that is perpen- 
dicular to the bow shock surface. Then our ex- 
pression for dEM is given by 



dEM = dS M 



(18) 



The differential surface area for any function w — 
f{z) used for a surface of revolution about the z- 
axis is given by 



dS — 271 w 



1 



duj 
dz 



1/2 



dz 



(19) 



Recalling that g = dz/dw =, T/Ta = {l+g'^)'^, 
and defining t = T/Ta, the rate change in sur- 
face area with the monotonic temperature "coor- 
dinate" along the surface becomes 



dS_ 



iV/' dz 



'^-^' + ^) dT 



27rii7 



^/l+g^ dz dw 



-1/2 271 



g dm dT 
du7^ 



dT 



T ' (jrVD ). 

dr 



(20) 



This rather elegant result requires yet one more 
step to fully determine the surface gradient en- 
tirely in terms of temperature; that step requires 
the solution to w{T). The mapping between w 
and T is made through the factor g using equa- 
tions ([5]) and yielding the implicit relation 
of 



g = am 



1. 



(21) 



Solving for w, taking the gradient with tem- 
perature, and doing some algebraic manipulation 
yields finally 



dS_ 
d^ 



(1-r) 



-(m-2)/(m-l) 



(22) 



We know that the result of the simulation for 
the emission measure is dEM/dT oc j^jjjj 
we know that the value of m for the bow shock 
is only a little steeper than a parabola. Using 
equations (fT8|) and ([22]) . it is convenient therefore 
to express the final result for the emission measure 
gradient as 



11 



dEM_ ^ TTN^R^'Ai _ 

dT (am)2/("-i) Ta ^ ' ' 

(23) 

where a = {m — 2)/(m — 1). If the bow shock 
were exactly a parabola, then m — 2, a — 0, and 
dEM/dT (X is a pure power law with an 

exponent remarkably close to the value of —7/3 = 
-2.33. With TO = 2.34 and a = 0.25, the OTSh 
approximation predicts an emission measure gra- 
dient of dEM/dT OC T"2.25 _ t)-0-25^ ^j^g 

limit that T ^ Ta, we have that r <C 1 and 
dEM/dT is again a power law with an exponent 
that is now even closer to the fit value. Obtaining 
an exact match to the fit exponent of the simu- 
lation requires m = 2.20 in the low t limit, not 
far from the to = 2.34 value derived from the 
bow shock shape. It would appear then that the 
OTSh is an excellent approximation for the emis- 
sion measure and temperature distributions of the 
bow shock in the hypersonic and adiabatic limits. 

There is one problem with equation (|23[) in that 
near the bow shock apex, r approaches unity in 
which case dEM/dT becomes singular for a > 0. 
This singularity is however only a mathematical 
artifact. The cusp arises from our fit to the ^(ro) 
formula with a power-law exponent of 2.34, but in 
fact the analytic result of Lomax & Inoye indicates 
that very near the apex the curve is a parabola 
in which case dEM/dT is non-singular. Earlier 
we noted that the p — —4/3 law underestimates 
the EM very close to the apex temperature. It 
is interesting that the parameterization of equa- 
tion ([23]) does lead to a rise in EM{T) when t 
is not negligible. In practice the singularity can 
be avoided by applying a high temperature cut- 
off, in which case the OTSh provides an excellent 
prescription for z{vu), T{tu), and EM(T). 

With regards to the A£ factor, based on 
the simulations, achieving the good results for 
dEM/dT in terms of S{T) means that the range 
of constant density along the bow shock is consis- 
tent with A£ a constant width. We find a value 
of A£ « 0.1i?c, independent of the y location 
along the bow shock. It is now useful to define an 
emission measure scaling parameter EMo, using 
a = 0.35 and m = 2.34, giving 



_ 7:N^R,'Ai 
° ^ (am)2/(™-i) 

where we assume a strong shock such that N = 
4:Nw, and Rc and N\y have both been scaled to 
the values used in the simulation. Using A£ = 
0.1i?c, the scale constant becomes EMo ~ 7 x 
lO'^o cm-3. As shown by WC07, typical EM val- 
ues for line formation are « 10^^ cm""^, implying 
that about 10^ clumps would be required. Assum- 
ing simple Poisson statistics, the expected vari- 
ability would then be about 1%, consistent with 
the low levels of X-ray variability from OB stars 
that has been previously noted. Interesting with 
EMo OC Rc^ , even small changes in the clump ra- 
dius - certainly a poorly known quantity - has 
a significant impact on the predicted numbers of 
the clumps and expected levels of variability sug- 
gesting that X-ray variability is a useful means 
of constraining clump properties (e.g., Oskinova 
et al. 2001). 

Although the giants and MS stars with lower 
mass-loss rates (smaller optical depths) show more 
dispersion in the range of Rfir, there is no evi- 
dence of any X-ray emission arising from below the 
associated optical depth unity radii. Furthermore, 
regardless of luminosity class, none of the stars 
show high energy ion stages forming far from the 
star. In the clump bow shock picture this means 
that the relative speed associated with the shocks 
is well below the ambient wind speed, which would 
produce very high temperatures. Clumps far from 
the star are likely being dragged out by the wind, 
and this reduces the shock temperatures in the 
outer regions, but with a sufficient shock jump in 
speed to produce the relatively low ion stages such 
as O VII, and N vii. 

As for the observed near-zero velocity shift and 
broad HWHM of the lines, these can best be ex- 
plained if we have the X-ray line emission aris- 
ing from the sides of the star as seen by the 
observer (i.e., perpendicular to the line-of-sight). 
This region contains a spread in speeds around 
the, Vz — 0, iso-velocity surface. Such a concen- 
tration of emission from that sector of the wind 
would occur if there were "self absorption" occur- 
ring. That is, each cool clump tends to absorb 
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the hne emission originating at the bow shock at 
the star- ward side of the clump. Thus we see an 
enhanced contribution from the regions at, and 
symmetric about, = 0. This counters the ex- 
pected tendency to see the blueward-shifted and 
-skewed line emission. 

6. Summary 

One of the advantages of the clump bow shock 
picture that it allows for a consideration of a 3- 
dimcnsional model for the X-ray production from 
OB star winds. The bow shocks can be distributed 
randomly about the star or in some other pre- 
specified way. Similarly one could choose to as- 
sume a variety of clump sizes at any radial shell. 
So a variety of ideas can be tested with the bow 
shock picture. Accounting for bow shocks around 
clumps in winds can potentially explain a number 
of properties of the X-ray emission observed from 
hot stars. 

The wide range of ionization stages that occurs 
at essentially every radial distance from the star 
can be produced by bow shocks because of the 
power law distribution of EM{T). Low ion stages 
such as O VII occur far from the star, but not the 
highest ions. One might think that the shocks 
would be strongest there because the wind speed 
is the largest. However, clumps could be dragged 
out by the wind and thus there could be a lower 
relative speed at large r, as discussed by Howk 
et al. (2000). 

The zero centroid shift problem, {Vpeak ~ 0.0) 
can be explained if the line radiation is gener- 
ated primarily from clumps that are to the side 
of the star. This could occur because of the "self- 
absorption" effect, whereby the clumps along the 
line-of-sight to the star strongly attenuate the X- 
rays that are produced on the starward side of the 
clumps. 

The broad line widths or HWHM of up to about 
half the wind speed can be caused by a combination 
of the sideways velocity the bow shock flow plus 
a range in values for an ensemble of clumps. 
This is in contrast with planar shocks that have 
been proposed for the shock fragments in hot star 
winds. The transverse velocity in the bow shocks 
is an essential aspect of the adiabatic or expansion 
cooling of the bow shock gas. 



The narrower line of main sequence and gi- 
ant stars can be explained by the correlation of 
wind speed with the transverse velocity around a 
clump. These stars have optically thin winds and 
the emission can be dominated by shocks located 
deeper in the wind where the inflow speeds are 
slower. 

The lines of high ion stages are not narrow as 
one might expect from source regions close to the 
star. However, to a large part this can be ex- 
plained by the relatively low spectral (or velocity) 
resolution. For the Chandra HETGS/MEG spec- 
tra the HWHM resolution is about 270 kms~^at 
25 A for lines of low ion stages, and to 1300 km s~^ 
for the high ions. The transverse velocity behind 
a clump and the range in clump velocities would 
lead to further broadening. 

There are observed EM{T) power laws (e.g., 
WS05), and these result directly from the wind- 
oblique shock interactions. The shape of the 
EM(T) distribution depends on the fraction of the 
wind material lying above optical depth unity. 

The Rfir radius corresponds rather well to the 
radius of optical depth unity, and this can probably 
be explained by having the observed line formation 
sources located to the sides of the star relative to 
the line-of-sight to the observer and having the 
mass loss rate reduced by an acceptable amount 
such as 1/3 the traditional M values. 

In addition to a consideration of the problems 
summarized above, we have introduced the OTSh 
picture. This provides a way to derive the line 
source information and the line of sight Doppler 
velocity shifts needed to compute line profiles from 
a single clump or an ensemble of clumps dis- 
tributed through a wind. 

The clump bow shock picture forms a useful 
element in a more realistic 3-D view of the stel- 
lar winds of hot stars, and the OTSh approach 
forms a useful technique for predicting and analyz- 
ing X-ray emission properties. The results should 
be of interest to areas other than hot star astron- 
omy. Clump flow interactions are thought to oc- 
cur in such objects as RS Oph (Nelson et al. 2008). 
Chimps with bow shock interfaces with a wind are 
seen in Hubble images of planetary nebulae (Odell 
et al. 1995). An understanding of clumps can also 
affect research well beyond the stellar wind com- 
munity. 



13 



We thank Nick Mocckcl for early work on the 
EM{T) relation and Greg Tracy for producing 
some of the figures. We are also grateful to John 
Brown for many interesting conversations on the 
topic of wind clumps. JPG, RI, AL and NAM have 
been supported in part by award TM3-4001 issued 
by the Chandra X-ray Observatory Center. We 
also acknowledge support from the NSF Center for 
Magnetic Self Organization in Laboratory and As- 
trophysical Plasmas. WLW acknowledges support 
by award GO2-3027A issued by the Chandra X- 
ray Observatory Center. Chandra is operated by 
the Smithsonian Astrophysical Observatory under 
NASA contract NAS8-03060. 

REFERENCES 

Abbott, D. G., Bieging, J. H., & Churchwell, E. 

1981, ApJ, 250, 645 

Barlow, M. J. and Cohen, M. 1977, ApJ, 213, 737 

Bouret J.-C, Lanz T., Hillier D. J. 2005, A&A, 
438, 301 

Bouret J.-C, Lanz T., Hillier D. J., Heap, S. R., 
Hubeny, I., Lennon, D. J., et al.2003, ApJ, 595, 
1182 

Brown, J. C, Ignace, R. & Cassinelli, J. P. 2000, 
A&A,356, 619 

Brown, J. C, Teller, D., Li, Q., Hanuschik, R., 
Cassinelli, J. P.,& Kholtygin, A. 2004, MNRAS, 
252, 1061 

Burke, A., Ignace, R., & Cassinelli, J. P. 2006, 
BAAS 38, 1046 

Canto, J. & Raga, A. 1998, MNRAS 297, 383 

Cassinelli, J. P. & Abbott, D. C. 1980, "The Uni- 
verse at UV Wavelengths NASA pub. 2171, 
pl27 

Cassinelli, J. P. Castor, J. I., & & Lamers, 
H. J. G. L. M. 1979, PASP, 90, 496. 

Cassinelli, J. P., Cohen, D. H., MacFarlane, 
J. J., Drew, J. E., Lynas Gray, A. E., Hoare, 
M. G.,Valerga, J. V., Welsh, B. Y., Vedder, 
P. W., Hubeny, I., & Lanz, T. 1995, ApJ438, 
932 

Cassinelli, J. P. Cohen, D. H., MacFarlane, J. J., 
Drew, J. E., Lynas Gray, A. E., Hubeny, I., 
Valerga, J. V., Welsh, B. Y., & Hoare, M. G. 
1996, ApJ460, 949 



CassineUi, J. P. & Hartmann, L. . 1977, ApJ, 212, 

488 

CassineUi, J. P. & Olson, G. L. 1979, ApJ, 229, 
304. 

CassineUi, J. P. & Swank, J. H. 1983, ApJ, 271, 
681. 

Cho, J. & Lazarian, A. 2002, Phys. Rev. Lett. 
88, 245001 

Cohen, D. H., Leutenegger, M. A., Grizzard, 
K. T., Reed, C. L., Kramer, R. H., & Owocki, 
S. P. 2006, MNRAS, 368, 1905 

Cox, A. N., ed., 2000, Allen's Astrophysical Quan- 
tities (4th ed.; Springer- Verlag) 

Donati, J. -F., Howarth, I. D., Jardine, M. M., 
Petit, P., Gatala, G., Landstreet, J. D., Bouret, 
J. -C, Alecian, E., Barnes, J. R., Forveille, T., 
Paletou, F., & Manset, N. 2006, MNRAS, 370, 

629 

Donati, J.-F., Babel, J., Harries, T. J., Howarth 
L D, Petit, P. & Semel M. 2002 MNRAS 333, 

55 

Evans, C. J. Crowthcr, P. A., FuUerton, A. W., 
Hillier, D. J., 2004, ApJ, 610, 1021 

Feldmeier, A., Kudritzki, R.-P., Palsa, R., Paul- 
drach, A. W. A., & Puis, J. 1997b, A&A, 320, 

899 

Feldmeier, A., Puis, J., & Pauldrach, A. W. A. 
1997a, A&A, 322, 878 

Feldmeier, A., Oskinova, L., & Hamann, W. -R. 

2003, A&A, 403, 217 

FuUerton, A. W., Massa, D. L., & Prinja, R. K. 
2006, ApJ, 637, 1039 

Groenewegen, M. A. T., Lamers, H. J. G. L. M., 
& Pauldrach,A. W. A. 1989, A&A, 221, 78 

Harnden, F. R., Jr., Branduardi, G., Gorenstein, 
P., Grindlay, J., Rosner, R., Topka, K., Elvis, 
M., Pye, J. P., & Vaiana, G. S. 1979, ApJ, 234, 

L51 

Hayes, W. D. & Probstein, R. F. 1966, Hyper- 
sonic Flow Theory, Volume I, Inviscid Flows. 
Academic Press. 

Hamann, W.-R., Feldmeier, A. & Oskinova, L. M., 
2007 Clumping in Hot Star Winds Potsdam 
University Workshop proceedings, in press. 



14 



Henrichs H., Neiner, C, & Geers, V. C. 2003 lAU 
Symp. 212, A Massive Star Odyssey, from MS 
to SN{ed. K. van der Hucht) ASP Conference 
Series, 305, 301. 

HiUier, D. J. 1991, A&A, 247, 455 

Hirschi, R. 2007, to appear in Potsdam Workshop 
Proceedings, in press (astro-ph 0709.0392) 

Ignace, R., & Gayley, K. 2002, ApJ, 568, 594 

Howk, J. C., Cassinelli, J. P., Bjorkman, J. E., & 
Lamers, H. J. G. L. M. 2000, ApJ, 534, 348 

Lamers, H. J. G. L. M. & Morton, D. C. 1976, 
ApJS, 32, 715 

Lamers, H. J. G. L. M. & Rogerson J. 1978, A&A 
66 417 

Lepine, S. & Moffat, A. 1999, ApJ 514,909 

Li, Q., Cassinelli, J. P., Ignace, R., Miller, 
N. A., & Brown, J. C. (2008) ApJ in press, 



( [astro-ph/0710263p ) 

Liu, X. & Osher, S. 1998, J. Comp. Phys. 141, 
561 

Leutenegger, M. A., Paerels, F. B. S., Kahn, S. 
M., & Cohen, D. H. 2006, ApJ, 650, 1096 

Lomax, H. & Inouye, M. 1964, Numerical Analysis 
of Flow properties about a Blunt Body NASA 
tech. rept, R-204 

Lucy, L. B. 1982, ApJ, 255, 286 

Lucy, L. B., & Solomon, P. M. 1970, ApJ, 159, 879 

Lucy, L. B., & White, R. L. 1980, ApJ, 241, 300 

Lupie, O. L., & Nordsieck, K. H., 1987, Astron. J, 
93, 214 

MacFarlane, J. J., CassineUi, J. P., Welsh, B. Y., 
Vedder, P. W., Vallerga, J. V., & Waldron, 
W. L. 1991, ApJ, 380, 564 

Macfarlane, J. J. & CassineUi, J. P. 1989, ApJ, 
347, 1090. 

Massa, D., FuUerton, A. W., Sonneborn, G., 
Hutchings, J. B. 2003, ApJ, 586, 996 

Miller, N. A., CassineUi, J. P., Waldron, W. L., 
MacFarlane, J. J., & Cohen, D. H. 2002, ApJ, 
577, 951 

Moeckel, N., Cho, J., & CassineUi, J. P. 2002, 
BAAS 34, 772 



Nelson, T., Orio, M., CassineUi, J. P., StiU, M., 
Liebowitz, E., & MucciareUi, P. 2008, ApJ, in 
press (astro-ph/0709.4512) 

Nugis, T. & Lamers, H. J. G. L. M. 2000, A&A, 
360, 227, 279 

Odegard, N., & CassineUi, J. P. 1982 ApJ, 256, 
568 

Oskinova, L. M., Feldmeier, A., & Hamann, W.-R. 
2006, MNRAS, 372, 313 

Oskinova, L. M., Feldmeier, A., & Hamann, W.-R. 
2004, A&A, 422, 675 

Oskinova, L. M., Ignace, R., Brown, J. C, 
CassineUi, J. P. 2001, A&A, 373, 1009 

Owocki, S. P., Castor, J. 1., & Rybicki, G. B. 1988, 
ApJ, 335, 914 

Owocki, S. P., & Cohen, D. H. 2006 ApJ, 648, 565 

Seward, F. D., Forman, W. R., Giacconi, R., Grif- 
fiths, R. E., Harnden, F. R., Jr., Jones, C, & 
Pye, J. P. 1979, ApJ, 234, L55 

Smith, R. K., Brickhouse, N. S., Liedahl, D. A., & 
Raymond, J. C. 2001 ApJ, 556, L91 

ud-Doula, A., & Owocki, S. P. 2002, ApJ, 576, 413 

Waldron, W. L. & CassineUi, J. P. 2001, ApJ, 548, 
L45. 

Waldron, W. L. & CassineUi, J. P. 2007, ApJ, ApJ, 
668, 456 (WC07) 

Wojdowski, P., & Schulz, J.P. 2005, ApJ, 616, 630 
(WS04) 

Wojdowski, P., & Schulz, J.P. 2005, ApJ, 627, 953 
(WS05) 



This 2-column preprint was prepared with the AAS lAT^rjX 
macros v5.2. 



15 



Table 1 

bow shock simulation parameters. 



Simulation Parameter 


Value(s) 


Rc 


10^" cm 


{Mach#)w 


47 & 71 


Vw 


1000 & 1500 km s"^ 


Nw,e 


IQio cm-3 


Tw 


2 X lO'^ K 


Ts 


1.4 X 10^& 3.2 X 10^ K 


PwVw 


100 & 233 dyne cm'^ 



Table 2 

Emission Measure Slopes versus Temperature 



Star 


Spectral Type 


Slope p 






{=l + d\ogEM/dlogT) 


CPup 


04f 


-1.1 


C Ori 


09.7 lb 


-1.2 


L Ori 


09 III 


-1.3 


S Ori 


09.511 


-0.9 


HD206267 


06.5 V 


-2.2 


P Cru 


B0.5 III 


-2.1 


T ScO 


B0.2 V 


-0.6 


Ori C 


04-6p 


+1.5 
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(a) 



(b) 



(c) 



(d) 




Fig. 1. — Simulation results of an adiabatic 
bow shock forming in response to an impenetrable 
sphere using the parameters given in Tab. [1] with 
Vw = 1500 kms~^. Shown are (a) flow stream- 
lines, (b) density with contours at 0.4, 1.2, 2.0, 
2.8, and 3.6Nw, (c) temperature with contours 
at 1, 3, 6, 10, 15, 21, and 28 MK, (d) emission 
measure contours of separated by factors 

of two with the largest value at 2.44 x 10^^ cm~^ 
and showing that the emission measure is isolated 
strongly at the shock front. 




z/io'°cm 

Fig. 2. — A contour plot of temperature in the 
bow shock simulation. The peak temperature con- 
tour value near the stagnation region is Ta given 
by eq. (121), and the contours in temperature have 
the same spacing as in Fig. [TJ Arrows show the 
velocity vectors. 



53.0 



52.5 



52.0 



51 .5 




6.0 



6.5 7.0 

log T (K) 



7.5 



Fig. 3. — The temperature distribution of the 
emission measure from a simulation of an adia- 
batic bow shock using the parameters given in 
Tab. [H The maximum temperature Ta is set 
by the speed of the incident wind relative to the 
shock. The low temperature cut-off is due to the 
finite size of the simulation grid. The dashed line 
is a power law approximation with a slope of —4/3. 
Two results are shown: these have Mw = 47 and 
71, which correspond to the I4ei= 1000 kms~^and 
Vrei=1500 kms~^as indicated 
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Fig. 4. — An illustration of the changes in velocity 
components at the shock front. The angles Ai 
and A2 are the incident and emergent angles of 
the flow relative to the direction perpendicular to 
the bow shock, and 6 represents the angle between 
the post-shock velocity vector and the z direction. 
The impact parameter of a streamline is given by 
its m coordinate value. 




ra/10'°cm 

Fig. 5. — The maximum temperature as a func- 
tion of impact parameter in our simulation of an 
adiabatic bow shock. Shown are simulation results 
(solid) and eq. [16] (dotted). 



